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Abstract. 

The most recent lattice results on the QCD phase diagram in the presence of a strong 

magnetic field strongly disagree with results from previous lattice simulations and several 

model calculations. The most remarkable difference is the qualitative behavior of the critical 

temperature as a function of the magnetic field intensity, which in later results are shown 

^ |. to decrease in opposition to what was found previously. According to the authors, such a 

O ' discrepancy in lattice simulations could be due to different lattice spacing, or different number 

i^H^, of flavors and quark masses. We investigate the influence of quark masses on the Polyakov- 

Quark-Meson model and show that, although quantitatively the results are sensitive to those 
parameters, the qualitative behavior remains the same. 

in 

^. Introduction 

^f-N , The phase structure of the QCD phase diagram in the presence of magnetic field was studied 

in recent years by several model calculations and also by full lattice QCD simulations. Such 
an interest on the subject was motivated by the fact that it was proposed that a very intense 
magnetic field, of order 5 — lOm^ [H [2], must be produced in non-central heavy-ion collisions 
and must have a considerable influence in the phase transitions suffered by the matter generated 
C^ . after the collision. 

Most of the model calculations O IU O [6] indicate an increasing behavior of the critical 
temperature as a function of the magnetic field and early lattice simulations ^ corroborated 
these results. Later lattice results [8], however, indicate a decrease of the critical temperature 
with the magnetic field intensity. In order to check their code and algorithm simulation, the 
authors of the latter have reproduced the results of the former group applying the same setup. 
Since in this context the previous results were recovered, the discrepancy was then associated 
to different lattice spacing or difference in number of flavors and quark masses. Approaches 
beyond mean field were considered [9], still obtaining an increase of the critical temperature 
with the magnetic field. A later approach [lUj could reproduce qualitatively the lattice results 
in [8], considering the context of the MIT bag model. 

According to [8] the difference between their result and the previous one can be traced back 
to the behavior of the chiral condensate as a function of the magnetic field. While in [7| it was 
found that the condensate has a monotonically increasing dependence on the field for all the 



range of temperature tested, [8] obtains a more complex behavior, in which such a dependence 
can increase or decrease depending on the temperature. In order to get some insight on the role 
of quark masses in model calculations we use the linear sigma model coupled to quarks and to 
the Polyakov loop to obtain a phase diagram as a function of the magnetic field for different 
values of constituent quark masses. Although we include the Polyakov loop, taking into account 
the interplay between the chiral and deconfining transition, we limit our analysis to the chiral 
condensate, following the indication in [8j that the chiral transition guides the general behavior 
of both transitions. 

The Polyakov- Quark- Meson model 

In this context, the deconfinement properties are encoded in the expectation value of the 
Polyakov loop, which works as an order parameter for this transition: 

l/T 

Confinement : s )t\ /r,' i ■ t m , Lix) = -TrPexp i / dr AAx^t 

[ (L) 7^ , high T ' ^ ' 3 L 7 



where j44 is the matrix-valued temporal component of the Euclidean gauge field A^ and the 
symbol V denotes path ordering. On the other hand, the expectation value of the chiral 
condensate works as the order parameter for the chiral transition, indicating if this symmetry 
is broken or restored for each value of external parameters: 

rhiral svmmetrv ■ f (^) / ' 1°^ T (/. = (a, vr) , 

Chiral symmetry . | ^^^ ^ q ^ j^-g^ ^ , ^ ^ (^+^^0^^-), (^) 

The chiral sector coupled to quarks is given by the Lagrangian: 

C = ^^[i^^^^^-g{a + i^^T■^)]^lJf + ]^{^^a^^'a + ^^^^^'^)-V{a,^), (3) 

which contains spontaneous and explicit symmetry breaking and the coupled of mesons and 
quarks via a Yukawa term. The quarks act as a heat bath in which the condensate evolves and 
quark loop corrections provides the thermal contribution to the potential. The implementation 
of the magnetic field is done replacing the simple derivative that describes the dynamics of the 
fermions by covariant derivative. Such a derivative contains an Abelian gauge field which can 
be seen as a vector potential for the magnetic field. Also the interaction with the pure gauge 
sector, through the Polyakov loop, can be implement in the same fashion, the Polyakov loop 
being parametrized in such a way that a non-Abelian gauge field interacting with the mesons 
can represent it. The 1-loop quark correction is then given by: 



e 



iV-M^JT 



det(ii^(g) - mq) 
det{i(^ — mq) 



det(ii/K9) _ ruq) 



(4) 



where the quark mass is given by rUq = \/m^ + [g{(j))'^ , with rric the constituent mass, and the 
covariant derivative contains both gauge fields, the Abelian and the non-Abelian. 

To reproduce the nature of the phase transition predicted in [SJ, a crossover, we take into 
account the quark degrees of freedom in the vacuum. The non inclusion of this correction implies 
a first order phase transition and a discussion regarding this point is done in [1]. 

We consider several values for the constituent quark mass in a range between the chiral limit 
and the physical consituent quark mass and also an extrapolation of that, rUq = 500 MeV, in 



order to investigate the systematics. For a set of temperature values, we obtain the condensate 
as a function of the magnetic field. The results are plot in Fig.l . In the chiral limit and for 
small values of the quark mass the condensate increases with the magnetic field and saturates 
at a constant value for all values of the temperature. This behavior is again in accordance with 
[7], which shows in their Fig. 3 a similar saturation (regarding that the periodic behavior is an 
artefact from the fact that they have introduced the U{1) field as a phase factor Ae^^). As the 
quark mass is increased the value of magnetic field to reach saturation becomes larger, however 
the behavior is systematically the same. 
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Figure 1. Expectation value of the chiral condensate for different values of magnetic fields. 
The different plots differ by the value of the quark masses considered and the different curves 



by temperature. We considered (a) ruq = MeV, (b) niq = 300 MeV and (c) 



m„ 



500 MeV. 



Conclusions 



In this work we have studied the role of quark masses for the chiral phase transition in the 
context of the Polyakov-Quark-Meson model in the presence of a magnetic field. Motivated by 
discrepant lattice results regarding the behavior of the chiral condensate as a function of the 



magnitude of the magnetic field we have studied the effect of going beyond the chiral hmit, 
since this was presented as a potential cause for the different lattice results. We showed that 
for all values of quark masses considered, from chiral limit to physical constituent quark masses, 
the qualitative behavior doesn't change, in accordance with [7J and in disagreement with more 
refined results in |8] . A systematic saturation of the value of the chiral condensate for all values 
of temperature is found, also in agreement with [7]. As a consequence, the critical temperature 
increases with the magnetic field and it was seen no sign in favor of a decreasing for any value 
of the magnetic field. 
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